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Abstract. The true nature of the radio emitting material observed to be moving relativistically in quasars and microquasars 
is still unclear. The microquasar community usually interprets them as distinct clouds of plasma, while the extragalactic com- 
munity prefers a shock wave model. Here we show that the synchrotron variability pattern of the microquasar GRS 1915+105 
observed on 15 May 1997 can be reproduced by the standard shock model for extragalactic jets, which describes well the 
long-term behaviour of the quasar 3C 273. This strengthens the analogy between the two classes of objects and suggests that 
the physics of relativistic jets is independent of the mass of the black hole. The model parameters we derive for GRS 1915+105 
correspond to a rather dissipative jet flow, which is only mildly relativistic with a speed of 0.60 c. We can also estimate that the 
shock waves form in the jet at a distance of about 1 AU from the black hole. 
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1. Introduction 

Microquasars are high-energy binary systems exhibiting jets 
with apparent superluminal motion suggesting that they are 
the miniature replicates in o ur Galaxy of the distant quasars 
; jMirabel & Rodriguezll99"8l GRS 1915+105 was the first mi- 
, croauasar discovered JMirabel & Rodrigueziri994l) and can be 
considered as an archetypical object of its class. Its radio 
emission can be divided into three distinct states: a plateau 
jet, pre- and post -plateau flares and radio oscillation events 
, jKlein-Wolt et alJ 120021) . The flares are the most powerful 
events. Their emission arises from two emitting blobs ob- 
served to move relativistically in oppos ite directions up to 
about 3000 AU from the core jFender et alJ 1 19991) . This 
has been interpreted as evidence for pairs of plasma clouds 
symmetrically ejected on b oth sides of the accretion disc 
(Rod riguez & Mirabellll999l) . This plasmoid ejection model 
is also used to explain the less powerful oscillation events 
arising in a shorter jet extending only up to ~ 30-50 AU 
JPhawan et aljEoOOl) . This interpretation is supported by the 
observation of simultaneous X-ray dips, suggesting that mate- 
rial from the disappearin g inner accretion disc is ejected from 
the black hole vicinity JBelloni et alJll997albl iMirabel etail 
ll998HEikenberrv et alJl998l) . 

Observational evidence for this phenomen on has rec ently 
also been obtained for the active galaxy 3C 120 jMarscher et alJ 



120021) . This discovery, by strengthening the analogy between 
galactic and extragalactic objects, rises the question on the true 
nature of the synchrotron emitting material. Is the microquasar 
plasmoid mode l or the quasar shock wave model the correct 
interpretation ? lAtovan & Aharonianl J 19991) showed that ex- 
panding clouds of relativistic plasma cannot describe the ob- 
servations unless there is an in situ electron acceleration pro- 
cess. An alternative model of a variable continuous jet also fails 
to explain the observed flux oscillations in GRS 1915+105 as 
it lea ds to unphysically large electron densities JCollins et alJ 
120031) . There is therefore growing evidence that the emitting 
electrons must be reaccele rated al ong the jet, as it is the case in 
a shock wave mechanism (Kai ser et all2 000). 

To further test the shock wave interpretation, we apply here 
to GRS 1915+105 the same model used to describe the long- 
term variability b ehaviour of 3C 273 from submillimeter to ra- 
dio wavelengths ( TiMex^LSlllQ QOj). T his model derived from 
the model of Marscher & Gear (1985) describes analytically 
the evolution of synchrotron emission resulting from a shock 
w ave moving down a r elativistic jet. Compared to the model 
of Ivan der Laanl J 19661) for an expanding plasma cloud it has 
the advantage to include the effects of both synchrotron and 
inverse-Compton energy losses. 

2. Data and Model 
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The data we use here are the radio oscillat ion events of 
GRS 1915+105 measured on 15 May 1997 JMirabel et alJ 
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Fig. 1. Model fit to the infrared and ra- 
dio observations of GRS 1915+105 on 
15 May 1997. The points in the four 
panels show the time dependence of 
the emission at observed wavelengths 
of 2.2yum, 2cm, 3.6cm and 6cm, from 
top to bottom. The 2.2 /im data are 
dereddened using a K-band absorptio n 
of 2.2 mag JChapuis & Corbell Eo04). 
The continuous red line is the best fit 
model, it is the sum of the emissions 
from the approaching jet (blue line) and 
the receding jet (green line) plus a con- 
stant flux at 2.2 fim fixed at 21.8 mJy. 
Red dotted lines separate the emission 
from distinct events. In the top panel, 
the onset times of the outbursts are in- 
dicated by arrows and the grey line is 
the X-ray light curve in the 2-60 keV 
band in units of 500 count s _1 . 



Il998t bhawan et alJl2000l) . The Very Large Array (VLA) ob- 
servations have the advantage to be very well sampled at three 
different radio wavelengths (see FigQ. The data start with a 
flux decrease followed by successive outbursts spaced by typi- 
cally 1-2 hours. These oscillation events are interpreted as ejec- 
tions in a small jet that was observed simulta neously by the 
Very Long Baseline Array (VLBA) (see iDhawan et alJ l2000. 
Fig. 9). Their evolution from short to long wavelengths is rem- 
iniscent of the observed behaviour of the b right quasar 3C 273 
on a time scale of years (Tiirl er et all2000l) . 

Additional data from the United Kingdom Infrared 
Telescope (UKIRT) in the 2.2 /mi infrared band show a flux 
enhancement which is very like ly the synch rotron precursor 
of the last radio outburst (M irabel et alll998l) . The simultane- 
ous X-ray measurements by the Rossi X-ray Timing Explorer 
(RXTE) in the 2-60 keV band present two flares apparently as- 
sociated with two consecutive radio outbursts. This X-ray be- 
haviour is indeed quite different from the light curve observed 
on 9 September 1997 which dis plays very rapid o scillations 
followed by a clear X-ray dip JMirabel et alJ fl998). A possi- 
ble interpretation of this difference is that the X-ray flares of 
15 May 1997 are not emitted by the accretion disc, but arise in 
the jet because of inverse-Compton scattering of synchrotron 
photons. 

The iMarscher & Geail ill 9851) shock model assumes a 
steady relativistic jet flow that can be disturbed by any increase 
of the pressure of the jet plasma or of its bulk velocity. If this 
disturbance is related to the accretion disc, it is likely to be sym- 



metric in both jets. Starting as a sound wave, the disturbance 
will become supersonic at some point because of the decreas- 
ing pressure along the jet. The resulting shock wave will prop- 
agate down the jet and accelerate particles crossing the shock 
front. Accelerated electrons will then emit synchrotron radia- 
tion in the increased magnetic field of the compressed plasma 
behind the shock front. The evolution of the emitted spectrum 
depends on the dominant energy loss mechanism of the elec- 
trons. Just after the onset of the shock, inverse-Compton radi- 
ation is likely to be dominant, then synchrotron losses should 
become more important and finally adiabatic expansion, due 
to the widening of the jet, will dominate the energy losses. As 
a consequence, the turnover of the self-absorbed synchrotron 
spectrum will move from high to low frequencies following a 
characteristic three-stage path (see Fig.|5Jl. 

3. Method 

To apply the model outlined above to the observed variability 
of GRS 1915+105 we use a meth odology similar to the one 
used previously for 3C 273 dTiirler et alJ2000l) . There are how- 
ever several modifications mainly implied by the specificity 
of GRS 1915+105, but including also some more general im- 
provements. An important difference between the modelling of 
3C273 and GRS 1915+105 is that for t he latter the emission 
from the receding jet is well observed dMirabel & Rodrigued 
ll994tlFender et al.ll999t lRodrfguez & Mirabell 19991) and thus 
cannot be neglected. We take this into account by modelling 
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Fig. 2. Evolution of the synchrotron emission of an average 
model outburst in GRS 1915+105 on 15 May 1997. The ob- 
served evolution is different for the approaching jet (a, c) and 
the receding jet (b, d). The grey lines are synchrotron spectra at 
times spaced by 0.2 dex in a, b, and equal flux density contours 
spaced by 0.3 dex in c, d. The thick black line is the path fol- 
lowed by the turnover of the spectrum. The triangles show the 
position of the two stage transitions for each outburst as com- 
pared to their average position (circles). The frequency of the 
four light curves in Fig.^are indicated by vertical lines. 

each outburst as being the sum of twin outbursts, one arising in 
the approaching jet and the other in the receding jet. The ob- 
served emission from the two opposite jets is however different 
due to orientation and relativistic effects. We model this by two 
parameters, one being the ratio between the Doppler factors of 
the emitting material in the two jets and the second being the 
observed delay between the onset of the twin outbursts. In addi- 
tion, the angle of 66 ° between the jet axis and the line of sight 
jFenderet all 1 999) is such that the shock is likely to be viewed 
sideways f or the ap proaching jet and from behind for the re- 
ceding jet dMarscher et al.lll992h . This was taken into account 
in the modelling and results in a different observed evolution of 
the two jets (see Fig.|5Jl. 

To define the specificity of each outburst, we allow to vary 
from one event to th e other only two ou t of the three parameters 
we used for 3C273 jTiirler et alJl2000|) . These two parameters 
are the normalization K on of the electron energy distribution 
(N{E) = K E~ s ) and the magnetic field strength B on at the on- 
set of the shock. The bulk Doppler factor £) on of the emitting 
material is assumed here to remain constant. Another simplifi- 
cation is that there does not seem to be a constant radio emis- 
sion component in GRS 1915+105 similar to the hot spot at the 
far end of 3C 273's jet. Furthermore, the start of the radio light 
curves of GRS 1915+105 can be simply modelled by the final 
decay of an outburst with average properties. 



To have a more realistic model we replaced by progres- 
sive transitions the low- and high-frequency breaks of the syn- 
chrotron spectrum (see iTiirler et"al] 12000. Fig. 2) and we 
smoothed out the sharp edges of the three-stage evolution. The 
normalization of the synchrotron spectrum was redefined as 
the extrapolation of the optically thin slope, as in the original 
iMarscher & GeaJ dl985h model. Finally, we introduced a new 
parameter to allow the frequency ratio Vh/v m of the low fre- 
quency break to the turnover of the spectrum to vary with time. 
This allows to describe a synchrotron source which is first in- 
homogeneous and becomes progressively homogeneous near 
the transition to the final adiabatic expansion stage, as it seems 
to be the case in GRS 1915+105 (see Fig.|2i,b). 

The model presented here has a total of 34 free parame- 
ters. 13 parameters are used to describe completely the shape 
and the evolution of the synchrotron spectrum for an average 
outburst in the approaching and the receding jets. The remain- 
ing 21 parameters are used to define the start times and the 
specificity of each distinct event. This model is the simplest 
jet model reproducing well the observations. It assumes a con- 
stant Doppler factor along the jet and a linear decrease B oc R l 
of the magnetic field B with the jet half-width R, as expecte d 
if B is perpendicular to the jet axis (Begelm an et alJfl984l) . 
The model parameters are constrained by a simultaneous least- 
square fit to the 795 radio and infrared measurements of Fig.^ 
This is achieved by many iterative fits of small subsets of the 
34 parameters. The overall fit is shown in Fig.^ Although the 
adjustment is not perfect (reduced^ 2 of 6.4), the shape of the 
light curves is very well reproduced by the model. 

4. Results and Discussion 

The overall evolution of an average outburst in GRS 1915+105 
as derived from the observations of 15 May 1997 is shown in 
Fig Et, c for the approaching jet and in Fig.|2j>, d for the reced- 
ing jet. The outburst is assumed to be intrinsically the same in 
both jets, but the observed evolution is different because of ori- 
entation effects. The different viewing angle changes the nor- 
malization, as well as the slopes of the three-stage evolution. 
The first effect is due to a different relativistic Doppler factor, 
while the second is due to the assumption that the approach- 
ing shock wave is viewed sideways, while the receding one is 
viewed from behind. 

The model parameters describing the jet properties suggest 
that the jet when observed was rather dissi pative and this coul d 
explain why it did not propagate very far (Dhawan et al. 2000). 
A first indication for this is the non-linear increase of the radius 
R oc L r of the jet opening with distance L along the jet. The 
best fit value of r — 1 .52 accounts for the rapid change of the 
spectral turnover frequency with time (see Fig. |2j;, d). It sug- 
gests that the inner jet of GRS 1915+105 is more like a trumpet 
than a cone. This shape could result from a pressure gradient 
of the material surrounding the jet or from diverging external 
magnetic field lines channeling the jet plasma. Alternatively, 
it might be the signature of a decelerating jet flow. To test this 
hypothesis, we relax the constraint of a constant Doppler factor 
along the jet. This results in a model with a slight tendency to- 
wards deceleration, but the improvement of the fit is not signif- 
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icant due to the lack of further infrared or submillimeter obser- 
vations. There are however other indications suggesting a dis- 
sipative jet flow. In particular the fact that the obtained value 
of s — 2.02 for the index of the electron energy distribution 
N(E) - K E s is close to 2. This means that there is roughly 
as much total energy carried by high-energy electrons than by 
low-energy electrons. For s < 2, the strong radiative losses of 
the dominant high-energy electrons would lead to a substantial 
pressure decre ase along the jet a nd prevent the shock to propa- 
gate far ( Marscher & G earll985l) . That such radiative losses are 
not negligible is also suggested by the slightly steeper decrease 
(k = 3.10) of the normalization K oc R~ k of the electron energy 
distribution than expected if the jet flow was really adiabatic 
(k = 2(s + 2)/3 = 2.68). 

From the ratio of the relativistic Doppler factors in the re- 
ceding and the approaching jet, we can derive the bulk speed 
of the emitting region. We obtain a ratio of 0.61, which corre- 
sponds, by assumi ng an angle of 66 ° between the jet axis and 
the line of sight dFender et all 19991) . to a speed of 0.60 c, c be- 
ing the speed of light. This mildly relativistic value contrasts 
with the speeds exceeding 0.9 c derived for the g iant eruptions 
dFender etalJll999t iRodrfguez & Mirabell ll999). but is com- 
patible with the hypothesis of a shock viewed sideways in the 
approaching j et provided that the j et opening half-angle is of 
13 ° at least JlVlarscher et ai1ll992l) . We can also estimate the 
distance from the black hole at which the shock wave forms 
based on the time interval between the onset of the outburst in 
the two jets. This interval is found to be ~ 400 s and corre- 
sponds to a distance from the black hole measured along the 
jet of 1.47 10 13 cm, by assuming the same orientation as above. 
This is roughly 1 AU and is only abou t twice the separation of 
the binary system jGreiner et al.l200ll) . 

Concerning the specificity of individual events, represented 
by the points in Fig. [2] we note that they do mostly differ in in- 
tensity from the average outburst. According to the model, this 
is due to changes of the normalization ^T on of the electron en- 
ergy distribution at the onset of the shock. The effect of signifi- 
cant changes of the magnetic field strength B on from one event 
to the other would lead to a clear scatter of the points along 
the frequency axis, which is not observed (see Fig. |2^, b). The 
fact that the specificity of individual outbursts is mainly due 
to differences in the normalization K on of the electron energy 
distribution is to be expected if the outbursts are due to an in- 
creased injection rate at the base of the jet. This finding there- 
fore supports the idea that the shock waves forming in the jet 
of GRS 1915+105 are the result of instabilities in the inner ac- 
cretion disc. 

It is difficult to assess whether a scaled-up version of the 
shock model presented here could apply to the giant eruptions 
of GRS 1915+105. Their radio l ight curve with its very fast 
rise compared to its slow decay i Dhawan et alJl200ol) resem- 
bles the infrared model light curve for the approaching jet, but 
expanded about 1000 times in duration and 100 times in ampli- 
tude. The optically thin synchrotron spectrum observed in the 
radio during the decay of major eruptions has also roughly the 
same slope than the infrared spectrum deduced here. This anal- 
ogy suggests that giant eruptions could be qualitatively sim- 
ilar events than those studied here, but with physical condi- 



tions such that their emission is shifted by orders of magnitude 
towards lower frequencies, longer time scales and enhanced 
fluxes. 

5. Conclusion 

This work shows that the microquasar GRS 1915+105 has, at 
least during some epochs, a variability behaviour similar to the 
quasar 3C273, which can be modelled by shock waves propa- 
gating in a relativistic jet. This analogy suggests that the phys- 
ical nature of relativistic jets and their emission is not depen- 
dent on the mass of the black hole at the origin of the ejection. 
The study of jets in both quasars and microquasars is therefore 
complementary. While quasars are bright and their jet struc- 
ture is relatively well resolved, microquasars have the advan- 
tage of shorter variability time scales and of a less relativistic 
flow allowing us to observe also their counter-jet emission. The 
ability shown here to constrain jet emission models by multi- 
wavelength observations opens new perspectives for studying 
relativistic jets in quasars, microquasars and tentatively even in 
gamma-ray bursts. 
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